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Critical exponents of the three-dimensional classical plane rotator 
model on the sc lattice from a high temperature series analysis 
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High temperature series expansions of the spin-spin correlation function for the plane rotator (or 
XY) model on the sc lattice are extended by three terms through order (3^"^ . Tables of the expansion 
coefficients are reported for the correlation function spherical moments of order I — 0, 1,2. Our 
analysis of the series leads to fairly accurate estimates of the critical parameters. 



In three dimensions the two-component vector model is the simplest spin model in the universality class of the 
] superfluid A transition of ^He, and of the ferromagnetic transition of magnets with an easy magnetization plane [ 
[01]. No high temperature (HT) series studies of this model have appeared in the last two decades in spite of 
remarkable experimental measurements of the critical parameters in superfluid "^He and intense theoretical activity in 
■ Renormalization Group calculations and by direct Monte Carlo simulations. 
I ' In particular we should mention that the critical index v, which describes the leading singularity of the superfluid 
Oh, fraction in ^He near the superfluid transition temperature, has been measured with high precision in a long series 
^ ' of experiments by G. Ahlers and his collaborators [ Ql]. As stressed by Ahlers, the superfluid fraction is the most 
accurately known singular parameter at a critical point, and correspondingly i/ is the most accurately known critical 
index. The most recent experiments yield the value v = 0.6705 ± 0.0006. 

Unfortunately the critical exponent 7 cannot be measured in liquid ^He and, as far as magnetic systems are 
concerned, no precise measurements exist either for 7 or for ly. A review of static critical properties of "^He can be 



o pre 

' found in Ref.[ |3|] and a general discussion of the interpretation of the measurements on ''He in connection with the 
problem of confluent singularities is given in Ref.[ [Q]. 

The Hamiltonian of the three-dimensional plane rotator (or XY) model is 

iJH--^ ^ s(x)-s(x + e^). (1) 

Here s{x) is a two-component classical spin of unit length associated to the site with position vector x — niei -\-niei -|- 
^363 = (fT-i, ?^2, ^^3) of a 3-dimensional simple cubic lattice and ei, 62, 63 are the elementary lattice vectors. The sum 
over X extends to all lattice sites. 

It has been rigorously proved that the model exhibits a ferromagnetic phase transition [ [||]. 

We present here series which extend by three terms, to order i3^'' , the series of Ref. [ [||]. They have been computed 
by a FORTRAN code which iteratively solves the Schwinger-Dyson equations for the correlation functions [ [Q]. 
We have tabulated the HTE coefficients of the two-point correlation function 

C{x;i3)^<s(fS)-s{x)> (2) 

for all inequivalent sites x for which the expansion is non trivial to order 0^'^ . 

We have analyzed the series for the spherical moments of the correlation function rffi^iji) defined as follows: 
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m«(/3) = ^|a;rC(x;/3)=^a«/3^ 

X r—1 



(3) 



(here \x\ — yjn\ + n\-\- n\ ), ? > and the sum extends over all lattice sites. The zeroth order spherical moment 
m}-^\f}) is the (reduced) susceptibility and is also denoted by xiP)- 

In Table I, we report the HTE coefRcients of the spin-spin correlation function < s(0) • s{x) > with x — (1, 0, 0). 

In Tables II, III and IV we report the expansion coefficients for the moments m^''\(3) with 1 — 0,1 and 2. 

Our analysis of this 0(2) symmetric model parallels that of the corresponding series for the 0(0) symmetric self- 
avoiding walk (s.a.w.) model and the 0(1) symmetric Ising model on the s.c. lattice [ ||]]. Using both first-order 
and second-order differential approximants, we first analyse the susceptibility series. We find that if the degree 
of the inhomogeneous polynomial is too low (< 3) the approximants cannot adequately accommodate the analytic 
background term. On the other hand if the degree of the inhomogeneous polynomial is too large (> 8), there 
are insufficient series terms to adequately represent the singular part of the series. For intermediate values of the 
degree of the inhomogeneous polynomial however, the approximants are stable, allowing the unbiased estimates 
Pc — 0.45406 ± 0.00005, 7 — 1.315 ± 0.009 to be made. The unbiased estimates from second-order approximants 
were more erratic, giving /3c — 0.4541 ± 0.0001, 7 — 1.32 ± 0.01. Biasing the approximants at (3c = 0.45406 gave 
7 = 1.315 ± 0.005 and 7 = 1.316 ± 0.005 from first-order and second-order approximants respectively. The results 
remain essentially unchanged using either the Fisher- Au Yang/Hunter-Baker definition (no regular singular point at 
the origin) or the Guttmann/ Joyce definition of the DAs (a regular singular point at the origin), the difference being 
mostly in the dispersion of the data (in particular of the background), which is somewhat greater with the former 
definition. 

A similar analysis of the first moment, m*^^)(/3) is less satisfactory. Most of the approximants are defective, but 
the few that are not are centred around a slightly lower temperature, Pc = 0.4542 ± 0.0003, with exponent 7 + = 
2.00 ± 0.03. Biasing the approximants at j3c = 0.45406 gives mainly defective approximants, slowly decreasing in 
value, so that we can only estimate ^ -\- v < 2.00. The second correlation moment series mP'\l3) is somewhat better 
behaved, though, like the analogous s.a.w. and Ising series, unbiased approximants at first glance give a lower critical 
temperature than do the susceptibility series approximants, notably /3c = 0.4542 ± 0.0002, and 7 -f 2i/ = 2.69 ± 0.02. 
This behaviour of the series mP'\(3) was also noted for the Ising and s.a.w. model series [ |^]. It appears that longer 
series are needed for higher moments of the correlation function. Biasing the second-moment series at /3c — 0.45406 
gives ^-\-2v — 2.67, but this must be regarded as an upper bound as the sequence of estimates of 7-|-2z^ decreases with 
increasing numbers of terms - just as observed previously for the corresponding Ising series. While it is difficult to 
extrapolate this slowly declining sequence, the limit 2.661q q2 hkely sufficiently conservative to include the correct 
value. In reaching this conclusion we have not only extrapolated this sequence, but have studied the behaviour of 
analogous sequences for the Ising and s.a.w. model, where we also have exact results for the two-dimensional models 
to guide us. If this estimate is accepted, we find from our earlier estimate of 7 that v = 0.67 ± 0.01. 

We may also construct the series with coefficients Cr = m^"^ /m^^ and study its singularity at z = 1 which should 
have exponent 2i'+l. Then we get v — 0.68±0.01, with again a decreasing sequence of exponent estimates, suggesting 
that f is in fact a little lower. 

In conclusion our estimates of 7 are fairly precise and, as shown later, in good agreement with the Renormalization 
Group (RG) results. However our estimates of 1/ cannot yet compete either with the precision of the experimental 
data nor with the RG or Monte Carlo determinations, although they are perfectly compatible with both. This is 
probably due to the slow convergence of m*-^''(/3), and as noted above, has already been observed in the study of 
high order expansions for the SAW and Ising models [ [^]. Longer series are then required in order to make a more 
accurate analysis possible and in particulare to account properly for the confiuent singularities. 

Let us now briefly review previous high temperature series analyses, restricting our review to the sc lattice results. 

Bowers and Joyce [ [0], computed series to order /3^ and gave the following estimates: Pc = 0.4530 ± 0.0016, and 
7 = 1.312 ±0.006. 

In Ref. [ ||ll|] the series were extended to order /3^^. The estimated inverse critical temperature was /3c = 0.4539 ± 
0.0013 and the corresponding estimates for 7 and v were 7 — 1.32 ± 0.05 and v = 0.675 ± 0.015. A comparison with 
our results shows that our central values for /3c and 7 are significantly lower and that the precision in our estimates 
has improved by a factor two. 

Reliable Monte Carlo simulations with good statistical accuracy, on reasonably sized lattices, have become possible 
only recently after the invention of algorithms with reduced critical slowing down [ |12|] . The largest accurately 
studied lattice is still only 64"^ sites large (present practical limits seem to be around 100^ sites), which means that 
a very accurate treatment of finite size effects is required and that the estimate of systematic errors is very delicate. 
The oldest analysis is due to Li and Teitel [ [|l^] who performed a Metropolis simulation ( supplemented by over- 
relaxation method ) on lattices up to 16^^ sites. A finite size scaling analysis of their data yields /3c — 0.4533 ± 0.0006 
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and V = 0.67 ± 0.02. (The model actually simulated is a clock model with 512 states.) 

More recently Hasenbusch and Meyer [ jlj]] used the WolfF single cluster algorithm on lattices up to 96^ sites. From a 
fit of the data to X oc (/3c-/3)"'^, they found /3c = 0.45421±0.00008 and 7 = 1.327±0.008. A recent update [ [|l5|] of this 
study using the WolfF single cluster algorithm on lattices up to 64'^ sites gave /3c — 0.45420 ±0.00002, v = 0.664 ±0.006 
and 7= 1.324 ±0.001. 

W. Janke [ jl^] also used the Wolff single cluster algorithm on lattices up to 48"^ sites. From a study of the fourth 
order cumulant he obtained /3c = 0.4542±0.0001 and v = 0.670±0.002. Fitting data to the formula x oc x+(/3c-/3)"^ 
he obtained /3c = 0.45408 ± 0.00008, and 7 = 1.316 ± 0.005. Repeating his fit with fixed /3c = 0.4542 the value of 7 
increases to 7 = 1.323 ± 0.002. 

The previous computations should also be compared to the estimates by the Renormalization Group applied to an 
0(2) symmetric 0'* field theory model. 

Sixth order perturbation expansion in three dimensions by Baker, Nickel and Mciron [ ||l^], gave 7 = 1.316 ±0.009 
and V = 0.669±0.003. Subsequently, taking into account the large order behavior of the perturbation series coefficients, 
Le Guillou and Zinn Justin [ |po[l] r efined these estimates and obtained 7 = 1.316 ± 0.0025 and ly = 0.6695 ± 0.001. 

Performing the computation [ |^l|] by the Wilson-Fisher e = 4 — d expansion Borel resummed to order e^, Le Guillou 
and Zinn Justin subsequently obtained the following estimates: 7 = 1.315 ± 0.007 and i' — 0.671 ± 0.005 . 

It thus appears that the RG results for 7 are slightly smaller than the old HT and some of the new Monte Carlo 
estimates, but perfectly compatible with the results of our analysis, while our estimate of is compatible with, but 
less accurate than, the most recent RG results. 



ACKNOWLEDGMENTS 



Our work has been partially supported by MURST. AJG wishes to acknowledge support from the Australian 
Research Council. 



[1] V. G. Vaks, A. I. Larkin, Sov. Phys. JEPT 22 , 678 (1965) 

[2] K. H. Muller, F.Pobell, G.Ahlers, Phys. Rev. Lett. 34, 513 (1975); D.S. Greywall, G.Ahlers, Phys. Rev. A7, 2145 (1973); 

L.S. Goldner, G. Ahlers Phys. Rev. B45, 13129 (1992). 
[3] P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 435, (1977), P. C. Hohenberg, Physica 109 , 110, 1436 (1982); 

H. Kleinert, Gauge fields m condensed matter, Vol. I, (Singapore, World Scientific 1989). 
[4] G.Ahlers, Rev. Mod. Phys. 52, 489, (1980), G.Ahlers, Phase transitions Cargese 1980, Edited by M. Levy, J. C. Le Guillou 

and J. Zinn- Justin New York, Plenum Press 1982. 
[5] J. Froehlich , T. Spencer Comm. Math. Phys. 83, 411 (1982) 

[6] M. Luscher and P. Weisz, Nucl. Phys. [FS22] B300, 325 (1988); P. Butera, M. Comi and G. Marchesini, Phys. Rev. B41, 
11494 (1990). 

[7] P. Butera, M. Comi and G. Marchesini, Phys. Rev. B33, 4725 (1986); P. Butera, R. Cabassi, M. Comi and G. Marchesini, 
Comp. Phys. Comm. 44, 143 (1987); P. Butera, M. Comi and G. Marchesini, Nucl. Phys. B300, 1 (1988); P. Butera, M. 
Comi and G. Marchesini, Phys. Rev. B40, 534 (1989). 

[8] D. S. Gaunt, A.J. Guttmann in Phase Transitions and Critical Phenomena, C. Domb and M. S. Green Eds., Vol. 3 
, (London, Academic Press 1974); G. A. Baker and D. L. Hunter, Phys. Rev. B7, 3346, 3377(1973); A. J. Guttmann in 
Phase Transitions and Critical Phenomena, C. Domb and J. Lebowitz Eds., Vol.13, ( London, Academic Press 1989); G. 
A. Baker, Quantitative Theory of Critical Phenomena, London, Academic Press 1990). D. L. Hunter and G. A. Baker, 
Phys. Rev. B19, 3808 (1979); M. E. Fisher and H. Au-Yang, J. Phys. A12, 1677 (1979); J.J. Rehr , G, S. Joyce and A.J. 
Guttmann, J. Phys. A13, 1587 (1980); B. G. Nickel and J. J. Rehr, J. Stat. Phys. 61, 1 (1990). 

[9] A. J. Guttmann, J. Phys. A20, 1839, 1855, (1987) 
[10] R.G. Bowers and G.S. Joyce, Phys. Rev. Lett. 19,630 (1967). 
[11] M. Ferer, M. A. Moore and M. Wortis, Phys. Rev.BS , 5205 (1973). 

[12] J. Goodman, A. D. Sokal, Phys. Rev. Lett. 56, 1015 (1986); R. H. Swendsen, J. S. Wang, Phys. Rev. Lett. 58, 86 (1987); 

M. Creutz, Phys. Rev. D36, 515 (1987); U. WoIfT, Phys. Rev. Lett. 62, 361 (1989); J. Goodman, A. D. Sokal, Phys. Rev. 

D40, 2035 (1989); A. D. Sokal, Nucl. Phys. B (Proc. Suppl.) 20, 55 (1991). 
[13] Y.H. Li, S. Teitel, Phys. Rev. B40, 9122 (1989). 
[14] M. Hasenbusch, S. Meyer, Phys. Lett. B241, 238 (1990). 



3 



[15] A. P. Gottlob, M. Hasenbusch, S. Meyer ( in Proceedings of the Int. Symp. on Lattice Field Theory "Lattice 92" Amsterdam 
1992). 

[16] W. Janke , Phys.Lett. A148, 306 (1990). 
[17] U. Wolff, Nucl. Phys. B322, 759 (1989). 

[18] R. G. Edwards, J. Goodman, A. D. Sokal, Nucl. Phys. B354, 289 (1991). 
[19] G. A. Baker, B. G. Nickel , D.I. Meiron, Phys. Rcv.B17, 1365 (1978) 

[20] J. Le Guillou, J. Zinn Justin, Phys. Rev. Lett. 39, 95 (1977) and Phys. Rev. B21, 3976 (1980) ; J. Zinn Justin Quantum 

Field Theory and Critical Phenomena, Clarendon Press Oxford 1989. 
[21] J. Le Guillou, J. Zinn Justin, J. de Physique Lett. 46, L137 (1985), J. de Physique, 48, 19 (1987). 



TABLE L HTE coefficients of the nearest neighbor corre- 
lation C(0, x) with X = (1, 0, 0) 



order 


coefficient 


1 


0.50000000000000000000000000000000 


3 


0.43750000000000000000000000000000 


5 


1.01041666666666666666666666666667 


7 


2.49169921875000000000000000000000 


9 


7.48240559895833333333333333333333 


11 


24.7292479338469328703703703703704 


13 


86.7042412409706721230158730158730 


15 


317.800753506891941898083560681217 


17 


1205.06602454131493586174488161069 


TABLE IL HTE coefHcients of the susceptibility m^°\ 


order 


coefficient 





1.00000000000000000000000000000000 


1 


3.00000000000000000000000000000000 


2 


7.50000000000000000000000000000000 


3 


18.3750000000000000000000000000000 


4 


43.5000000000000000000000000000000 


5 


102.343750000000000000000000000000 


6 


237.054687500000000000000000000000 


7 


546.946289062500000000000000000000 


8 


1252.00488281250000000000000000000 


9 


2858.81752929687500000000000000000 


10 


6496.15140787760416666666666666666 


11 


14735.3746412489149305555555555555 


12 


33314.753774685329861 1 111111111111 


13 


75222.2566392081124441964285714286 


14 


169444.488235923222133091517857143 


15 


381306.311343971793613736591641865 


16 


856543.263379992410619422872230489 


17 


1922537.91945074856684251367029620 
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TABLE III. HTE coefficients of the first correlation mo- 
ment m^^\ 



order coefficient 

~0 0.00000000000000000000000000000000 

1 3.00000000000000000000000000000000 

2 11.4852813742385702928101323452582 

3 35.3919166429113710288472410676167 

4 100.391645797382835211177404733391 

5 270.169140885332810622174619548742 

6 703.928165009702962171567107355945 

7 1789.19653133764917963889959830865 

8 4468.32789180460469625866305929854 

9 1 1000.87266696858 1 17344288428576 16 

10 26788.0560947846126416923232831814 

11 64627.3429637161982839763298977200 

12 154749.818273239775925845196634614 

13 368132.797893714045088109726930470 

14 870977.871997489895140365839695762 

15 2050710.75491296809029207572988208 

16 4808405.28831745065018387682551317 

17 11232374.4966970585972846903939187 



TABLE IV. HTE coefficients of tlie second correlation mo- 
ment m^^K 

order coefficient 

~0 0.00000000000000000000000000000000 

1 3.00000000000000000000000000000000 

2 18.0000000000000000000000000000000 

3 72.3750000000000000000000000000000 

4 247.500000000000000000000000000000 

5 770.593750000000000000000000000000 

6 2261.34375000000000000000000000000 

7 6360.66503906250000000000000000000 

8 17343.7773437500000000000000000000 

9 46158.4210449218750000000000000000 

10 120515.319303385416666666666666667 

11 309746.425031873914930555555555556 

12 785831.296427408854166666666666667 

13 1971809.99205790928431919642857143 

14 4901417.59164962163047185019841270 

15 12084656.3170853394364553784567212 

16 29584235.7640201335230832377438823 

17 71970593.8709586784015817546900548 
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